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Abstract 



Starting from a static spherically symmetric solution of the Einstein's 
field equations in the second approximation in the perfect fluid scheme, 
^ ' in the exterior of the source, the corresponding metrics depend on a di- 

mensionless parameter a . Taking the Sun for the source, the influence of 
a on the classical Solar System tests of the general relativity is estimated, 
the only a-dependent one being the Shapiro radar-echo delay experiment. 
Performing such a test in given conditions within a precision better than 



10 2 , it is possible to obtain an experimental value for a . If a = 1, the 



o 

m . 

, Einstein's equations in the weak field approximation take the D'Alembert 

■-. ■ form, which attests the existence of gravitational waves 1 

O . 
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■ For the metric theories of gravitation, the retardation of light signal passing 

near massive bodies provides a principal test (Shapiro, 1964-1989, Weinberg, 
1972, Reasenberg et al., 1979). In the following, we shall perform a more ac- 
curate analysis of this effect, starting from its " anisotropy" , i.e. the direction 
dependence of the light velocity in the gravitational field of the Sun - which dis- 
tinguishes it from other classical Solar System tests. The origin of the anisotropy 
is to be found in the manner of gauging the gravitational potentials, the effect 
being evidenced by means of a certain dimensionless parameter a. In the second 
section, a statically symmetric solution of the Einstein's field equations in the 
second approximation for a spherical mass distribution is given and the equa- 
tions of motion for both massive bodies and photons are derived. In the third 
section, the influence of this metrics on the Solar System tests of General Rel- 
ativity is estimated. The fourth section concerns with the connection between 
the Shapiro test and the propagation of gravitational waves. 
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2 The static field of the Sun and the equations 
of motion 



2.1 The static metrics 

We start from the standard static spherically symmetric metrics (Zeldovich, 
Novikov, 1971): 



(dS) 2 = e v(r \cdtf - e x{r) {drf + e 



l dn 



(1) 



where v{r), A(r) and cr(r) arc arbitrary functions to be established from phys- 
ical arguments. By asking this diagonal metric to fulfill the Einstein's field 
equations: 



H,v = 0,1,2,3 



(2) 



in the perfect fluid scheme (Misner, Thornc, Wheeler, 1973), wc come to the 
following expressions for the nonzero components of £^„ and T M „: 

E oa = [(Aa + fa' 2 ) + A') - ±AV + - 

Eli = - [\o 12 + \v'o' + + v>) + + ^e A " CT 
£ 22 = rV" A [-i(Aa + ia' 2 ) - |(Ai/ + ±V 2 ) + + A') + \a'{\' - v>) + \v'\'] 

E 33 = (sm 2 9)-E 22 

Too = e» [pc 2 +(pj^-p)' 

Tn = pe x T 22 = r 2 pe a T33 = (sin 2 6»)T 22 

(3) 

where p is the invariant rest mass density and p is the invariant pressure of the 
source fluid; we denoted by prime the derivative with respect to r . In the linear 
approximation, the field equations become: 



Aa+l(a'-\') + ±(a-\) = -^p 
Acr + Az/-i(i/ + A') = 



(4) 



Assuming the source to be a sphere of radius R , we have: 

G 
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M = — / PE{r) ■ 4:nr 2 dr A$ = -47r- 
c Jo < 



2 PE 



(5) 



1 This work has been partially supported by the 1033C Romanian Space Agency contract 
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where pe is the total rest energy density of the source, Mq is the total rest mass 
and $ the Newtonian potential. In the linear approximation, pe = pc 2 . In the 
second order approximation, the gravitational self energy as well as the elastic 
energy coming from the equilibrium pressure are included. The solution of 4 is: 



"2-2 A 



a — arbitrary 



(6) 



(see Tonnelat, 1965 about the arbitrariness of u). 

Outside the mass distribution, an accurate solution of the field equations 
= is: 



(dS) 2 



7W 



{cdtf 



// 



f{r) 



f 2 {r){drY + f 2 (r)dQ} (7) 



where, due to the asymptotic flatness, f(r) is subjected to the constraint: 



lim 



f(r) 



1 



r^oo T 



(8) 



The solution 7 may be checked by direct standard calculations. The integration 
constant p is determined as: 

GM 

V=— , (9) 

from the necessity of regaining Newtonian theory as a non-relativistic limit. 
From 8 we can write f(r) as an expansion: 



f(r) =r + ap + 0(^ r ) 



(10) 



provided that r and p arc the only arguments of / with physical dimensions. 
Here a is a dimensionless parameter. Using 10 we get an explicit approximate 
form of the metrics 7: 



(dS) 2 = 



1 - 2^ + 2a (cdt) 2 - [(l + 2^) (dr) 2 + (l + 2a£) r 2 dQ 



(11) 

A little discussion is here needed. Outside the source, derived from 1 and 
11, we have a w 2a^. On the other side, the quantities u, A, a cannot contain 
due to difficulties concerning the continuity conditions on the surface of the 
spherical source. The only possibility in hand is therefore a w 2a$ , which 
leads to: 



(dS) 2 ^ (l-2* +2a^ (cdt) 2 - 
- { [l + 2a$ - 2(1 - a)?f] (dr) 2 + (l + 2a$) r 2 dtt} 



(12) 
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The going through the solution 12 ensures the compatibility between the inner 
and outer solutions within the given approximation. Actually, it was proved 
that the metrics 12 is a second order approximate solution of the field equations 
2 provided that the rest energy density and the equilibrium pressure are: 

PE= (l-^^pc 2 + (a-l)($ + 2r$> + 2p p = - J p$' dr (13) 
Details are given in Appendix A. 



2.2 The equations of motion 

The equations of motion are resulting from the geodetic variational principle, 
which is equivalent to the Euler-Lagrange variational principle: 



i J Ldt = 







L = -m a c 2 K K = 



dS 
cdt 



v 



A-lB-r+C 



(v • r) : 



err 1 ' 



d (dL\ 



dt 



dL 
dr 



Here we used the energy integral: 

E = v • p 



dt 1 



dr 



K 



1/2 

(14) 
(15) 

(16) 



in order to eliminate the quantity ra^K^ 1 in the expressions of the momentum 
p and the force F , because moK ~ 1 becomes undetermined, 0/0, in the case of 
the photon. From 14 and 15, we get: 



P = 



1 raoc 2 d 

2 K dw 



K 2 



1 m c 2 d 

2 K dr 



(K 2 ) 



(17) 



As dE/dt = 0, the equation of motion acquires a form which holds for both 
massive bodies and photons: 



d_ 

dt 



1 d (\ 



Adv \2 



-K' 



1 d (\ 



Adv\2 







(18) 



Comparing dS/cdt from 14 to the homologous quantity from 11, we get: 



(19) 



Ak l-2 t - + 2a(^) Bwl + 2a- C«-2(a-l)- 
r \ r / r r 

By using series expansions relying on jit/r < 1, we finally derived the accelera- 
tion experienced by a point-like body of rest mass m , m <C M n , outside the 
source: 

09 T) 2 " 



aw — 



GMq 



1-2Ci- + 2C*2/3 2 -2C 3 - 



r + 2C 4 ^(/3-r)/3 (20) 
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where: 

d=a + l C 2 = l~ C 3 = -|(a-l) C 4 = a + 1 0=^(21) 

3 The Solar System general relativistic tests 

From the results of section II, performing usual computations (Weinberg, 1972), 
we derived formulas for estimating the general relativistic tests, within the Solar 
System, for an arbitrary gauge condition. The results are shown in the sequel. 

The perihelion advance of Mercury per revolution 5u> may be estimated by a 
standard procedure. The equation 20 is integrated in order to obtain the energy 
and angular momentum integrals: 



(22) 



( i v 2 _ GMi) + ^[(2C 4 - 2C 2 + I C 3 )v 2 + IC 3 ^ 
+(Ci + 2C 2 - 2C 3 - 2C 4 )]^ = £ 
(1 + 2C 4 ^) (vxr)=A 

£ = A- E -- U c 2 + , A=^M 

m y2 1 2 mgc 2 J ' E 

where is the constant areolar velocity, E the energy and M the angular 
momentum. Out of these equations, a relativistic Binet equation may be derived 
(Ionescu-Pallas, 1980): 

^ + k.( 1 -)=± + 2C^cos*6 (23) 
al9 2 \rj p! p 2 



where k = l+2-£-(Ci-2C 2 -2C 4 ), ft = GM /c 2 , Po = a(l-e 2 ), Pl = Po +O(fi). 
Then, a perturbational technique is used for reaching the trajectory formula: 

I = _L[ 1+ecos (^)] + | Ca ^( 1+sin ^) . (24, 
After such calculations, 5u) turns out to have the expression: 

5bJ = ( - Cl + 2 ° 2 + 2Ca) cMI-^) = c 2 a(l e 2 ) (25) 

which is independent on a. Here Mq is the rest mass of the Sun, a - semi major 
axis of the ecliptic of Mercury, e - Mercury ecliptic eccentricity. 
The light deflection angle 5^ is defined as: 

1 f + °° 1 
(5* = — / a x (t)dt= — M+oo) - ^(-oo)] (26) 
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where: 



a x = -^R\(l + 2C 2 )-2C 3 cos 2 9} r = (c 2 t 2 + R 2 V cos9 = %■ 
r R 

(27) 

Here, R is the Sun radius; a x is obtained from 20 by replacing there v ^cj and 
x = R, y = ct. v x is the projection of the photon velocity on the Ox axis (trans- 
verse to the trajectory) (Ionescu-Pallas, 1980). After standard calculations, we 

2 „ \ 2GM AGM 



5* = 1 + 2C 2 



c 2 R 



c 2 R 



that is (5 is also independent on a . 

Obviously, the gravitational red shift is independent on a , due to its con- 
nection only to the linear approximation of goo • 

However, the situation turns not to be the same with the Shapiro retarda- 
tion effect (see for the following Shapiro, 1964, 1989, Weinberg, 1972, Misner, 
Thornc, Wheeler, 1973). We shall derive the formula for ST - the retardation 
of a radar signal in a to-and-fro Terra-Mercury journey, when Mercury is at the 
upper conjunction. The elementary duration spent by a photon is, from 11: 



cdt 



1 + 2^ + (a 
r 



rdr 



Vr^ 



R 2 



+ R 



d( 9 - go) 
dr 



dr 



(29) 



The photon is travelling, let's say, almost along the Oy axis, passing at the 
nearest distance R from the center of the Sun. The first term in 29 leads to the 
retardation due to the gravitational refractive index n along the non-deviated 
trajectory (x = R) 

(30) 



n=l + 2^ + {a-l)R 2 ^ 



The second term leads to the "overtaking" retardation. There are no problems 
to estimate the first effect. For the estimation of the second effect, we first notice 
that 0o(r) = arccos j- is the equation of the line x — R (the trajectory in the 
absence of the gravitational field). Thereafter, we make use of the trajectory 
equation: 

1 COS 9 \l r 

in order to get the formula: 



(a + 1) cos 9 + (a - 1) cos 2 9] 



(31) 



(9 - 0„ ) 



R \pr 



R 2 



2-(a + l)- + (a-l)f- 
r \ r 



(32) 



whence 



R~r{9 — 9 ) = 
dr y U) Vr 2 - R 2 



a + 1 
r + R 



(a-1) 



R 



(33) 



G 



The integration is straightforward. We obtain for the retardation ST the for- 
mula: 



ST = 



4GM 



In 



4M 2 
R 2 



[220.5+ 19. 7{a + 1)] x 10~ 6 s (34) 



where l\ is the major semi- axis of the ecliptic of Terra and I2 is the major semi- 
axis of the ecliptic of Mercury. Performing a Shapiro-type experiment within 
an error of < 10~ 2 it is possible to determine an experimental value of a. 



4 The gauge parameter a and the gravitational 
waves propagation 

We write the general metric tensor as: 

= Viiv + V ( 35 ) 

where r}^ = 2<5 0/J <W — 5^ is the covariant Minkowski tensor. From 11 and 19, 
we have, in Cartesian coordinates: 

h Q0 = (A-l) hjk = -(B-l)S jk -C^ , j,k = 1,2,3 (36) 

Imposing the well-known Hilbert gauge condition (Zeldovich, Novikov, 1971, 
Tonnelat, 1965): 

1 \ 
h x °--rf a h),\ = A, a = 0,1, 2, 3 ( ), A = ) (37) 

where r} Xcr = 2<5° A <5 0cr — 5 Xa is the contravariant Minkowski tensor, with: 

h = i) Xa h x<J w 4o^ (38) 
r 

we get from 36, 37 and 38, by replacing the coefficients A, B, C with their values 
from 19: 

a = 1 (39) 

i.e. the condition for a in order to satisfy the Hilbert gauge. 

From the other point of view, Einstein's field equations 2, written in rectan- 
gular coordinates and in an inertia frame: 

R»v = -^- (t» v - ± 9llv T) (40) 

reduce, with T M „ in the perfect fluid scheme, and within the linear approximation 
35 (| VI < 1). t0: 



8ttG 



-P \ Uf,u v - -r^, 



(41) 
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where: 




(dS M ) 2 = rjxadx x dx a 




(42) 



1 a? 



For bringing 40 to the D'Alembcrt form, we need just the condition 37, i.e. a = 
1. The D'Alembcrt form unambiguously attests the existence of gravitational 
waves, in the framework of Einstein's General Relativity Theory. 

5 Conclusions 

We performed a more accurate analysis of the Shapiro time-delay Solar System 
test of general relativity. In order to perform the analysis, in the second section 
of the paper we studied a general spherically symmetric static metrics given by 
the Sun, in the second approximation. Starting from a more general solution 
of the Einstein's equations 7 and imposing the physical conditions of regaining 
the Newtonian limit and the continuity of the solution on the surface of the 
spherical mass source, we obtained an approximate form of the metrics in the 
exterior of the mass source 11, which contains, in the second approximation an 
undetermined dimensionless parameter a. By means of standard calculations 
performed using the metrics 11, we got a form of the equations of motion 18 for 
point masses, which holds for both massive bodies and photons, and we derived 
the acceleration 20 experienced by a point-like massive body or photon, both 
18 and 20 holding for the movement in the static field of the Sun. In the third 
section, we estimated the influence of the undetermined gauge parameter a on 
the classical Solar System tests. The perihelion advance, the light deflection 
and the gravitational red shift are independent on a. However, the Shapiro 
effect depends on a. Performing the calculations for the retardation of the 
radar signal in a to-and-fro Terra-Mercury journey, with Mercury at the upper 
conjunction, we estimated a relative influence of a of the order ~ 10~ 2 on the 
total signal retardation of time 34. Finally, in the fourth section, we established 
the connection of the parameter a to the propagation of the gravitational waves. 
From the Hilbert gauge condition 37, necessary for the Einstein's equations in 
order to take the D'Alembert form (which attests the existence of gravitational 
waves), we obtained the value a = 1. 

The short conclusions of our analysis arc the following: i) it is possible to 
determine an experimental value for a by means of a Shapiro-type experiment, 
and ii) if a = 1, the existence of general relativistic gravitational waves is proved. 

Finally, a definite alternative - which may be solved experimentally through 
the agency of the Shapiro effect - is reached: either the light propagation in 
the gravitational field of the Sun is isotropic and the gravitational waves do 
exist, or the respective propagation is anisotropic and the gravitational waves 
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are precluded. It is to point out that our conclusions are conditioned by the 
hypothesis about the correctness of Einstein's General Relativity Theory. 
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6 Appendix 

All the physical quantities entering the Einstein's field equations of General 
Relativity Theory are expanded in power series of the constant 1/c 2 : 



c 

X =^o + ±\i + O (1) (44) 

*=^o + ^i+oQf) (45) 

P = Pm- ^Pm Q a o + °o) +0^^ (46) 

P = Pm + o(-L) (47) 

(the label M stands for the Minkowskian nature of the respective physical quan- 
tity) . Accordingly, the tensors E a p and T aj 3 are themselves expanded in similar 
power series: 

^ = ?*3 + ?*S+0(^) ( 48 ) 

T Q ^c 2 r^+T«+0(l) (49) 

and the field equations E a p = — (8irG/c 4 )T a p split in a set of equations which 
no more contain the constant c: 

BQ = -**cn«$ , (50) 

i5g=-87rGT$ ,etc... (51) 

One obtains: 

T S = PM(r)6 Oa 5 Q0 , (52) 

T i/3 = \pm{t) \vq(R) + v {r) - A (r) - 2a (r)] 5 0q (5 0/3 - a af3 p M (r) (53) 

where a a( g is the tensor of Minkowskian metric in spherical coordinates, and 
Pm(i") is the equilibrium pressure: 

Pm{t) = p M {r)<S>' M {r)dr (54) 

J r 
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expressed by means of the Newtonian potential: 

A r $ M = -4ttGpm (55) 
The equations 50, written in explicit form, are the following: 

A r cr + ^ {cr'o - K) + ^ (°o - A ) = -8irGp M {r) 

-JK + ^o) + ^(Ao-<to)=0 (56) 

-^A r a - ^A r i/ + K + *o) = 

Assuming for er a value er = 2a$M, one obtains the solution of the system 
56 as: 



is Q = -2$ 



M 



A = 2a$ M (r) + 2(a - l)r$^ 



M ! 



a = 2a$ M (57) 



The system 51 contains not only the quantities u\,Xi,a\, but the quantities 
^o,Ao,cto as well. Inserting in 51 the values 57 of the zero-labelled quantities, 
one obtains 51 in explicit form, namely: 



12 

M 



[A r ai + i (a[ - X[) + £ (at - Ai)] + (a 2 + 2a - 2) $ 
8ttG {(a - 1) [$m - (a + l)r$' M ] p A f + [^Pm + Pm] } 



+ (a 2 - 2a + 2) = -87rGp 



[-|A r ai - |A r i^ + ^ K + A'J] - (a 2 - 2a + 2) $f f = 

8ttG I i (a - l) 2 p M ■ r$' M - p M } 
We multiply 60 by 2 and thereafter add it with 59 and 58; one obtains: 

A r ^+4 (1 - a) $ft = -8nG |(1 - a) [2r$^ - $ M ] PAf + 



(58) 
(59) 
(60) 



GMo 



R 



-Pm - 2p 



AI 



Now, let us define a new function $ by mean of the formula: 



Accordingly, 



$ = $m + 



$ 2 



(a - 1) $ 2 M - -v x 



1 - 2^ + 2a^j- = 1 + -^vq + 



1 fl 



c 4 V2 



(61) 



(62) 



(63) 
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From 61 and 62 we derive the Poisson type equation for <f>: 



A r $ = -^%p E , (64) 



p E = (l~ ^JjPj Pmc 2 + (a - 1) ($ M + 2r$' M ) p M + %Pm 



(65) 



12 



